Observation of hexatic liquid vortex matter in YBa2Cu207_j 
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An atomic beam probe is used to study the structure and dynamics of 

quantized supercurrent vortex lines in YBa2Cu307_5 at temperatures up to 

within 0.7 K below Tc- Here we report the direct observation of a vortex 

configuration with sample wide bond-orientational order but only short range 

translational correlation. The data imply the existence of an intermediate 

'hexatic' vortex line liquid phase. We find that the hexatic liquid is in thermal 

equilibrium over a narrow temperature range below Tc and is quenched into 

an immobile hexatic glass at low temperatures. 

A magnetic field enters type II superconductors in the form of discrete supercurrent 
vortex lines, each carrying a quantum of magnetic flux. The microscopic structure and 
dynamics of vortex matter is not only crucial for the application of superconductors - since 
it limits the capacity for carrying electrical currents - but also has a fundamental bearing 
on diverse topics such as the structure of liquid crystals |]^ and the phase diagram of high 
density DNA mesophases 0. In clean high-Tc superconductors there is evidence for a first 
order melting transition from a low temperature rigid and crystalline vortex lattice to a high 
temperature mobile and disordered vortex liquid |@,||,PJ3l- -^"^^ '^^^^ of YBa2Cu307_rf, 
there are indications that the melting transition changes character at critical points both at 
high (ifucp) dlB well as low magnetic fields (-fficp) |1T0|JTT|| . The interesting region of the 



vortex matter phase diagram just below Tc has previously proved inaccessible to vortex- 
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imaging-type studies and there are few detailed results on the microscopic structure of the 
vortex liquid. 

In analogy with theories for 2D melting [|12|], it is predicted that extended line objects such 
as nematic liquid crystals |jTB| and quantized flux lines |n| can exhibit a hexatic line liquid 
phase with long range bond-orientational order but short range translational correlation. 
The experimental determination of a hexatic phase hinges on two issues: (i) the comparative 
range of translational and bond-orientational order and (ii) whether the structure is in 
thermal equilibrium. Early observations of hexatic correlations in the vortex array by low 
temperature decoration techniques ||l5| , p!6| could not resolve whether the images represented 
a frozen-in vestige of a hexatic vortex liquid disordered due to thermal fluctuations or the 
ground state of a lattice disordered by static random pinning defects IT^, such as found in 



charge density wave systems [|T8|. Decoration studies of BSCCO samples indirectly indicated 



an equilibrium hexatic structure significantly below Tq [0. Recent Hall probe imaging of 
the vortex liquid phase in BSCCO pO| have been too limited to provide data on the range 
of correlations of the vortex liquid. Common to all these earlier reports is that they have 
probed a limited number of vortices (< 10^). We have used a new sensitive atomic beam 
technique to study directly vortex correlations in a high quality YBa2Cu307_5 single crystal 
in the low density regime below H\cp for temperatures to within 0.7 K below Tc- We have 
probed correlations in a sample containing 3 x 10^ vortices, and report evidence for a hexatic 
vortex line liquid with sample wide bond-orientational order. We further demonstrate that 
the hexatic liquid is in thermal equilibrium over a narrow temperature range below Tc and 
is quenched into an immobile hexatic glass at low temperatures. 

Experimentally, to probe vortex correlations, we pass a beam of potassium (^^K) atoms 
near the sample and detect magnetic resonance (MR) transitions driven by the microscopic 
magnetic field due to vortices (see Fig. la). Transitions are resonantly driven in atoms when 
the atomic velocity and the spacings of vortices along its path combine to make an oscillating 
magnetic field component at the magnetic resonance frequency /q. We actually measure 
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the excitation probability for atoms as a function of their velocity and work backwards to 
infer spatial characteristics of the vortex array. If the atomic beam passes over an ordered 
vortex lattice, the MR transition will be strongly driven for velocities v = dfo corresponding 
to vortex row spacings d along the beam (fig lb). Formally, to first order, the velocity 
dependent MR excitation probability, Pmr{v), represents a projection of the vortex matter 
structure factor PH : 



1 r 2 
Pmr{v) ~ ^ J dqy [h{ci,z,X)] S'2(q) , Qx = 27r/o/w. (1) 

Here h{q^,z, A) is the Fourier transform of the field from a single vortex with penetration 
depth A at height z above the surface ||22| and q = {q^, Qy) is the in-plane wave vector. In 
analogy with x-ray scattering, [h{q,z, A)]^ may be identified as the vortex form factor, and, 
most importantly, S'2(q) = / exp(— iq • r)g{r)d'^r is the 2D vortex lattice structure factor 



- the Fourier transform of the vortex pair distribution function g{r) ||23|. Details on the 
experimental method are published elsewhere pT| , p^ . 

The experiments we report here were performed on a 0.7 x 1.7 x 0.1 mm^ YBa2Cu307__5 
single crystal {Tc = 93.0 K, ATc < 0.3 K). The sample was grown by a self-fiux method 



|25| and thermo-mechanically detwinned. Two electrical contacts (< O.l^l) are placed 1.3 
mm apart on the top surface. The atomic beam skims the bottom surface along the a-axis 
and probes a 0.7 x 1.0 mm^ area. A typical MR scan takes 30 s to acquire. 

Figure 2 shows MR profiles recorded during cooling at a rate 0.2 K/min in a field 12 Oe 
parallel to the crystalline c-axis. We find that the MR profiles show two broad peaks for 
all temperatures T < 92.3 K for which we detect a signal. Initially, upon cooling, there is 
a gradual screening of the magnetic fiux from the sample. Below 90 K, the fiux density is 
constant at 10.7 G, corresponding to a lattice constant Oq — 1.5 fim. This screening is seen as 
the slight increase of MR peak distances when the sample is cooled. The two peaks represent 
(projected) vortex row spacings for a vortex 'lattice' with a unit cell that is misaligned with 
respect to the atomic beam (similar to Fig. lb, upper panel). In contrast, either an isotropic 
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vortex liquid or a 'powder pattern' of randomly oriented crystallites without sample wide 
orientational order would show a single bump representing the projection of a diffraction 
ring (Fig. lb, bottom panel) 0,0]. An immediate conclusion is that the vortex array has 
sample wide orientational correlation in the entire temperature range probed. 

A detailed analysis of the MR signal is necessary in order to extract the translational 
correlation length. The analysis is outlined here, details will be published elsewhere pi 



We model the atomic excitation using equation (1) and correct for saturation of the MR 
transition in a rate-equation approach [^]. Saturation is important since it broadens the 
MR peaks, especially when the magnetic field modulation is strong: close to the sample 
surface and at low temperatures (short penetration depth A). We model the vortex pair dis- 
tribution function g{r) = E«,m / (r - r^.^^, o" (n„)) , where f{r,(r) = {27Ta^y^ exp{-r^ /2cr^) 
and the sum runs over a lattice = ^ri + ^i'2 spanned by ri and r2 (except I = m = 0). 
Orientational disorder is included by averaging the structure factor 5*2 (q) for a range of 
primitive lattice vectors ri and r2 rotated over an angular interval of width A6. The com- 
puted transition probability (1) is corrected for saturation and averaged over the relevant 
range z < 2 fim for direct comparison to the experimental data. The main model freedom 
is in the functional form of the displacement 'dispersion' function a{r) and the orientational 
spread A6. The unit cell parameters ri and r2 are determined from the length scales for 
the two MR peaks, their relative intensities, and the flux density B, which is measured 
independently using a miniature Hall array. 

The best fit to the experimental line shapes (smooth curves in Fig. 2) is for a triangular 
lattice that is tilted 18° with respect to the atom beam direction and compressed along 



that direction — the crystalline a-axis of the sample p6|. The fit is for a sample wide 
orientational disorder A6 < 2°, a dispersion function of the form a{r) = cro(r/ao)^ with 
(To = (0.14 ± 0.01)ao and roughness exponent p = 0.5 =F 0.1. This dispersion implies a 
liquid-type short range order. The data could not be fitted using a quasi-long-range-order 
logarithmic dispersion predicted for the Bragg-glass phase |27]. The magnitude of the nearest 



neighbor dispersion is consistent with the Lindeman criterion (o"o/ao ~ 0.1 — 0.15) for a 
hquid near its melting point. Our experimental data therefore shows a vortex array with 
a liquid-like translational correlation length = (16 ± 8)ao (determined by the criterion 
o'{^t)/O'0 = 0.5) that is much shorter than the orientational correlation length ^ 500ao 
(sample size). 

Does this vortex pair distribution function represent a hexatic phase with true long range 



bond-orientational order or does it reflect finite size effects in a vortex liquid with slowly 



decaying, yet short range, orientational order? Toner |^ argues that point defects affect 
translational order more than orientational order, and predicts that the correlation lengths 
should scale as = ^t/^o in the absence of orientational coupling with the underlying 
crystal. In this scenario, we would expect < lO^aQ. We estimate a limit ^g > IO^Oq for 
the orientational correlation length by assuming that the angular disorder increases with 
distance as rapidly as the translational disorder, i.e. A9{r) ~ (r/ao)°'^. Here ^g is defined 
by A6{^g) = 30°. Clearly, this scaling is not satisfied - the orientational order is long ranged. 
The structural data in fig. 2 is therefore consistent with a hexatic vortex array with true long 
range orientational order. The absence of scaling is evidence that long range orientational 
order is induced by a coupling to the underlying crystal |^. A further indication of this 



coupling is a preferred direction for the vortex array with respect to the crystalline axes p9 
we find repeated field cooling runs always produced vortex arrays with the same orientation. 

Over which temperature range is the hexatic configuration in thermal equilibrium? The 
experiment shown in Fig. 3 addresses this issue. First, we observe that sample wide orienta- 
tional order is destroyed by applying a strong current to the field cooled vortex array at 82 
K. This is seen as the disappearance of the two distinct diffraction peaks for currents exceed- 
ing 220 mA (Fig. 3 a). (Note: 220 mA is almost an order of magnitude below the critical 
current for the sample at this temperature.) The broad MR line shape seen for large currents 
is consistent with a vortex array without sample wide orientational order. The stability of 
the disordered vortex array after the current is switched off (Fig. 3a, bottom curve) shows 



that the transport current has introduced plastic deformations - the vortex array is not in 
thermal equilibrium at 82 K. Second, the sample is subsequently warmed and we observe 
the annealing of the disordered vortex array above 91 K (Fig. 3b); Specifically, we see the 
reappearance of the two distinct diffraction peaks - indicating long range bond-orientational 
order. This shows that the hexatic structure is in a thermal equilibrium (liquid) state over 
a ~ 2 K interval below Tc- 

For temperatures below 91 K the structure of the hexatic liquid is quenched into a hexatic 
glass. This happens when the dynamics of the vortex liquid and hence the evolution toward 
crystallization is stopped as thermal fluctuation energies become smaller than potential wells 
due to pinning centers. This scenario is supported by several experimental observations: (i) 
The minimum pinning energy is comparable to thermal energies at the glass formation 
temperature. The threshold current (/min = 120 — 220 mA at 82 K) for the onset of change 
in the vortex array structure (Fig. 3a), provides an estimate for the minimum vortex pinning 
force, Fmin = (polrain/w ^ 5x 10^^^ N. Here w = 0.7 mm is the width of the sample. Assuming 
a quadratic pinning potential with a range comparable to the coherence length ~ 30 A), 
the corresponding minimum threshold pinning energy is f/min/^B — -^min^/2A;_B ~ lO^K. 
(a) The amount of disorder does not change measurably between 90 K and 10 K. The 
best fit model comparisons (blue curves in Fig. 2) at 10 K and 91.5 K are for correlation 
functions g{v) with identical amount of disorder. The effect of temperature is included in 
the model line shapes only through the penetration depth dependence of the vortex form 
factor [/i(q, 2;)]^ p2[ in equation (1). The two curves are for Aab = 140 nm and 400 nm at 
10 K and 91.5 K respectively. 



Theoretically, the hexatic flux liquid should be bounded by two phase transitions 12,14 



At high temperatures, it should transform to an isotropic liquid by the proliferation of 
unbound disclinations. At low temperatures, the hexatic liquid should undergo a phase 
transition to a well ordered state (crystalline or Bragg glass |^^) by freezing out unbound 
dislocations. Our experimental data do not show any of these transitions. Assuming that 
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the hexatic-to-isotropic transition follows the established high field first order melting line 
||7|.|6|JT0[] , a simple extrapolation predicts Tm — 0.9998Tc at 12 G - beyond the range of our 
method. On the low temperature side, we have shown that crystallization is preempted 
by the formation of a pinned hexatic vortex glass. The existence of a low density hexatic 
vortex liquid near Tq implies that the single first order melting transition for intermediate 
fields separate into two transitions in the region below the lower critical point i^icp. The 
robustness of iJicp with respect to oxygen vacancies and weak point defects due to electron 



irradiation |Tl|], further suggests that a low field hexatic liquid phase could be a generic 
feature of the magnetic phase diagram of YBa2Cu307_5 below Hicp. We believe that hexatic 
order is imposed on the vortex liquid due a coupling to the underlying electronic structure 
of the sample. The range of existence of hexatic phases in other high-Tc materials may 
therefore depend on the anisotropy of the crystal structure. 
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FIG. 1. (a) Experimental principle: The time dependent magnetic field seen by atoms 
moving near a sample with superconducting vortices drives magnetic resonance (MR) tran- 
sitions in atoms. If the field has a strong Fourier component at the MR frequency /o, then 
the transition 1^2 will be strongly driven. For vortices spaced a distance d apart, the MR 
transition will be strongly driven for atoms moving with velocities near v — dfo . (b) Model 
predictions for two different vortex lattices. In both figures, the grey arrow indicates the 
direction of the atomic beam. For an ordered vortex lattice, the transition probability of 
the atoms, , has maxima at velocities corresponding to lattice row spacings along the beam. 
The upper example shows a 10.5 G hexatic lattice rotated 15 degrees with respect to the 
atomic beam - similar to the case reported in this work. The two row spacings indicated, di 
and ^2, are readily identified with the two peaks in the data. Note: a triangular lattice with 
long range translational order would give a similar signal but with much taller and narrower 
peaks. For contrast, the lower example shows model predictions for 10.5 G isotropic liquid 
vortex structure. In reciprocal space, the signal corresponds to a projection of 2D vortex 
structure factor (shown in the insets) onto the beam direction. 

FIG. 2. Atomic beam scans showing hexatic vortex matter structure. The MR transition 
probability versus length scale is measured during field cooling (at 0.2 K/min) in an applied 
field 12 Oe||c. The blue curves arc best fit model predictions for identical structure factors 
with Xab— 140 and 400 nm at T=10 K and 91.5 K, repectively. The presence of two distinct 
MR peaks indicate sample wide orientational order. Note: the scale for each individual 
curve has been adjusted to fit on the graph. The signal strength at 10 K are 10 times larger 
than those at 91.5 K. 

FIG. 3. MR profiles showing destruction of orientational order by a transport current 
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and subsequent restoration by annealing, (a) Current dependence of MR profiles. The 
sample was initially field cooled to 82 K in a 12 Oe field (top curve). MR transition profiles 
are shown for increasing currents (middle three curves). The bottom curve shows the MR 
profile 1 hour after turning off the 300 mA current. The broad bump at high current is due 
to a vortex array without long range order, (b) Annealing of the disordered vortex array. 
The top curve shows a MR profile due to the current induced disordered vortex array (panel 
a) after the current was turned off. The gradual appearance of two distinct peaks when 
annealing above 90 K shows the re-appearance of long range bond-orientational order. 
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